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Abstract 

The sum of A'^ sufficiently strongly correlated random variables 
will not in general be Gaussian distributed in the limit — > oo. We 
revisit examples of sums x that have recently been put forward as 
instances of variables obeying a g-Gaussian law, that is, one of type 
est X [1 — (1 — q)x'^]^^^^~'^\ We show by explicit calculation that the 
probability distributions in the examples are actually analytically dif- 
ferent from g-Gaussians, in spite of numerically resembling them very 
closely. Although g-Gaussians exhibit many interesting properties, the 
examples investigated do not support the idea that they play a special 
role as limit distributions of correlated sums. 
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1 Introduction 



The central limit theorem says that the sum of N sufficiently weakly cor- 
related random variables, when properly scaled with and in the limit 
N ^ oo, has a Gaussian probability distribution. In many places in physics, 
however, strongly correlated variables occur. Causes of strong correlations 
may be, for example, the presence of long-range interaction or memory ef- 
fects. On the basis of a generalization of the statistical mechanical formalism 
[1] it has been suggested that in certain strongly correlated systems such 
many-variable sums are distributed according to a g-Gaussian law. This law 
is a generalization of the ordinary Gaussian and arises, in particular, when 
one maximizes the generalized entropy which is at the basis of the theory. 
Pursuant to this idea g-Gaussians have been advanced to describe velocity 
distributions in a Hamiltonian mean-field model (HMF) of classical rotators 
[2J, in turbulent Couette- Taylor flow |3j, and in cellular aggregates [4|; and 
were also proposed for the velocity distribution of galaxy clusters [5] and the 
temperature fluctuations in the cosmic microwave background [6] . In at least 
the case of the HMF, divergent interpretations of Monte Carlo results have 
led to vivid discussions in the literature (see e.g., [71181191 flO]). 

For a single variable x the g-Gaussian probability distribution Gq{x) is 
given by 



where Cg is a normalization constant and the scale factor b controls the width. 
For — oo < g < 1, which will be the range of g values of interest to us, Gq{x) 
is bell-shaped and has its support restricted to the interval [— Xm, x^n], where 

= (|E| ^) ^ • The denominator 3 — g in fll.ip is there merely to guarantee 
that for g — i> — oo the support of the g-Gaussian remains finite. The special 
case G_oo is the block function equal to ^b^ on [—b~2^b~2]. For g ^ 1 ex- 
pression (11. ip reduces to the ordinary Gaussian Gi{x) = (6/27r)^ exp(— 

Independently of their applicability in statistical mechanics, g-Gaussians 
have several interesting mathematical properties. For example, the multivari- 
ate g-Gaussian [for which in (II. ip is replaced with a symmetric bilinear 
form ^^fj^^x^A^^Xt,] has marginal distributions that are g'-Gaussians with 
g'^g(see|II]). 

We will be concerned specifically with the question of whether g-Gaussians 
can emerge as probability distributions of sums of correlated random vari- 
ables. Attempts in the literature to demonstrate how such g-Gaussian distri- 
butions may arise, have consisted of constructing and numerically analyzing 
simplified model systems. A successful model should not only explain the 
origin of the g-Gaussian, but also express the parameter g in terms of more 
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common concepts such as interactions and correlations. These model studies 
leave for later investigation the question of whether the properties preas- 
signed to the model may actually be realized in physical systems. 

In this note we will examine two model systems taken from the recent 
literature [121 US] . Both feature sums that were put forward as candidates 
for being g-Gaussian distributed. We show by explicit calculation that in 
both examples the probability distributions are actually analytically differ- 
ent from g-Gaussians, even though very closely resembling them numerically. 
Notwithstanding the interesting properties exhibited by g-Gaussians, the two 
examples studied therefore do not lend support to the idea that these func- 
tions play a role as limit distributions of correlated sums. 



2 First example 

2.1 A set of strongly correlated variables 

Thistleton et al. [I2] consider what is probably the simplest imaginable in- 
stance of a strongly correlated system. It consists of N identically distributed 
variables mi, M2, ■ ■ ■ ,un such that the probability distribution Pi of each in- 
dividual Uj is the block function 

Pi(n,) = l, (-i<«^. <i), (2.1) 

and Pi{uj) = elsewhere. The Uj have, furthermore, strong mean-field type 
correlations, as expressed by their covariance matrix 

(2.2) 

r (j k). 

where < r < 1 . The authors of reference [12] are then interested in finding 
the probability distribution Vn{U) of the scaled sum 

N 

U = N-'Y,u, (2.3) 

j=i 

and in obtaining its limit V{U) = \im]\r^oo'PN{U). It will turn out in this 
mean-field problem that in order for the limit — > 00 to exist, the appropri- 
ate scaling is the one with A^~^ given in (12.31) rather than the conventional 
one with N~^^^ applicable to independent variables. It is therefore clear that 
we cannot expect the limits r —>■ and ^ 00 to commute. 

Determining Vn{U) requires the knowledge of the full joint probability 
distribution P/v(mi, U2, ■ ■ ■ ,un) of the uj. This distribution is fixed implicitly 
by the way the Uj are constructed [12] , namely, from a Gaussian distribution 
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that one knows how to generate numerically. Let the correlated Gaussian 
variables zi, Z2, ■ ■ ■ ,zn be distributed according to 



N N 



P^^-^^{zuZ2,...,ZN)=[{2nfdetM]-Kxp(^-lJ2T.'^^ik'^)- (2-4) 

j=i k=i 

Then the zj have the covariance matrix (zjZk) = Mjk, for which will take 

M,,= { , ^ ^ (2.5) 

where p G (0, 1] is a parameter. It follows that the elements of the inverse 
matrix are 

f (j = k), 

[ -P {j ^ k), 

with 

^ l + (iV-2)p ^ P .2 71 

(i-p)[i + (iv-i)p]' ^ (i-p)[i + (iv-i)p]- 

The marginal distribution P^^^^'^^Zj) of derived from (12.41) is a Gaussian of 
unit variance. The variables uj are now constructed [12j from the Zj according 
to 

Uj = [ 'dxP^^'^'^ix) 
Jo 

= ^erf(^)' (j = l,...,iV). (2.8) 

Then Uj has the block distribution fl2.ll) and is correlated to the other Uk 
according to fl2.2p with an r that may be determined as a function of p. 

On the basis of numerical evidence Thistleton et al. [12j conjectured that 
the probability distribution 'P{U) of the scaled sum U of these variables was 
a g- Gaussian with an index q related to p according to 

g=i^. (2.9) 
1 -p 

This example is sufficiently simple to allow for an exact analytic treatment, 
which we will present now. 
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2.2 The limiting law ViU) 

In this subsection we replace the block function Pi considered in (12. ip by 
an arbitrary symmetric function Pi{u). Let there be N correlated random 
variables Uj, j = 1,2, ... ,N having an arbitrary symmetric marginal distri- 
bution Pi{uj) of unit variance and having correlations (ujUk) = r for j ^ k. 
Variables uj having these properties may be obtained from the zj by 



dxPi*^^"''(x) = / dxPi(x), {j = l,...,N), (2.10) 
Jo 

which generalizes (12.81) and which we will abbreviate as 

u, = h{z,), {j = l,...,N). (2.11) 
The expression for the law Vn{U) follows directly from its definition, 

/oo 
dzi... dzN6{U - A^-^ J2 ^(-^i)) 
-oo ■ 

xexp(-i J]5^z,M-iz,). (2.12) 

j k 

The quadratic form in the exponential in (I2.12p simplifies according to 

j k j j 

Upon substituting (I2.13P in (I2.12p and introducing integral representations 
for exp(i/5(X]j ZjY^ and for the delta function, we find that the integrations 
on the Zj factorize and we get 

fOO J,, /"OO J \ 



Vj,{U) = riV(27r)^detMl-^ r^e-^^'^' T— e- 
^ ^ i-oo 7-00 27r 

) . AT 

^^Q-\{a+l3)z^+u{NpY/-^z+i\N-^h{z)\ _ {2.14) 



We use that a+(3 = 1/(1— p) and shift from 2 to a new variable of integration 

(l-p)-h(z-u{l-p){Npy^). (2.15) 



w 



N/B 



X 



Completing the square in the exponential in (12.141) produces a factor exp ( 
(1 — p)z^^) which may be taken outside the integral and whose A^th power 
combines with exp(— |A^z/^) to yield exp(— ^7//^) where 7 = A^(l — p)/[l + 
(A^ — l)p]. The asymptotic A^ dependence of this coefficient, namely 7 ~ 
(1 — p)/p, is the result of a cancellation of terms of order A^. Using that 



5 



det M = [1 + (A^ — l)p](l — p)^ ^ and rearranging some of the prefactors we 
transform fl2.14p into 



X 



1 roo , lu 

i / . i„.2 I :\ ^r-l^,/^„ ^ n — 



'27r J 



du;e-^"''+*^^"'^«"'+'^)v^)) . (2.16) 



At this point we perform an asymptotic expansion in powers of ^ of the 
right hand side of equation (12.161) . This gives 

1 „ s 1 /"OO POO 1 \ 



ViU) = hm ^ / dz/e-^" / — e 



x(l + ^y^-^e~^-^/^((^ + .)v/w)+0(Ar-)) .(2.17) 
We now define 

k{u)=U - r ^e-'^""' h{{w + u)^l^p) (2.18) 
J~oo V 27r 

and simphfy (12.17^ further to 

1 — /*oo 

Let z/*(f/) be the unique solution of 

k{iy,) = 0, (2.20) 

so that we may write 5[k{h')) = 5(z/ — z/*)|A;'(z/*)|~^ Then, in the general case 
under consideration here, our final result for the limit function V{U) follows 
from flCTjl as 



2^' 



(2.21) 

in which /c is defined by (12.181) . z/* by (I2.20p . and we have set k' = dk/du. 



2.3 Special case of a block function 

The limit function V{U) given by equations (I2.2ip . (I2.18p . and (12.200 can be 
found fully explicitly for the special case of the block function distribution 
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(12. ip studied by Thistleton et al. [121. According to equations (I2.8P and 
(12. lip , we have in that case h{z) = |erf(z/A/2), which when substituted in 
(|2:T8|) yields 

k{u) = f/- ierf(KZ/), 

z/, = K-^eTr\2U), (2.22) 

where erf~^ denotes the inverse error function. After straightforward sub- 
stitution of (I2.22P in the general expression (I2.2ip we obtain 

V{U) = (^) ^ exp (- '-^[err\2U) f), {-\<U<\). 

(2.23) 

This is our end result for the distribution V{U) of the scaled sum U in the 
special case studied by Thistleton et al. [12] • It is not the g-Gaussian that 
the authors speculated it to be. 



2.4 Discussion of V{U) 

For U ±1 we have that erf~^(2f/) — > ±oo and therefore V{U) = 
outside the interval [— One may verify explicitly the normalization 

1 /2 

jJ^"^V{U)dU = 1. We consider now successively the limits p — > and 
p — >■ 1. The inverse error function has the expansion [T3| 

erf-\2U) =nm+ l{rrWf + U^Wf + l(vr^f/)^ + + • • • 

(2.24) 

For p <^ 1 it suffices to keep the first term in this expansion and substitu- 
tion in (I2.23P shows that in that limit 'P{U) approaches a Gaussian of vari- 
ance p/{4:7r). This Gaussian limit is nevertheless somewhat curious: whereas 
it is true that for p ^ the variables become independent, the usual 
Gaussian distribution associated with their sum would require a scaling with 
j\r-i/2 instead of the present one with A^~^. Hence we have here the explicit 
demonstration of the non-commutativity of the limits p —>■ and N —>■ oo. 

In the other limit, p — > 1, one sees directly that V{U) tends to unity on 
the interval [— |, |]. This second limit is heuristically obvious since for p = 1 
the random variables are fully correlated. For < p < 1 the function V{U) 
interpolates between a Gaussian and a block function. Since the g-Gaussian 
does the same for 1 > g > — oo, one may ask, for given p, whether there 
is a corresponding g that gives the best fit. This will be done in the next 
subsection. 
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Figure 1: Solid line: the exact limit function V{U) given by equation ()2.23p for 
the case P = jq- Dotted line: the g-Gaussian approximation to this function for 
g = — I (see text). The difference between the two curves is of the order of the 
thickness of the dotted line and just barely visible to the eye. 
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2.5 Best g-Gaussian fit to 



We try to find the best g-Gaussian fit to solution fl2.23p by adjusting the 
values of q and the width h in (11. ip . There is no unique way of defining a 'best 
fit'. We will base ourselves on comparing the small-argument expansions of 
V{U) and Gq{x). The one of Gq{x) is immediate, 



oo 

Gq{x) = b^Cq exp [ - ^ 



n-1 



(2.25) 



The expansion of V{U) is obtained by substituting (I2.24p in (12.231) and going 
through the algebra. One finds 



r{u) 



2-p\h 



p 
p 



exp 



oo 



n=l 



exp 



2(1 -P) 
P 



C2„7r U 



n=l 



(2.26) 



in which the coefficients are defined by the C2n which in turn are given 
by 

r — ^ r— 2 „_26 „_176 r — 8138 (0 07\ 

C2 — i, C4 — o, Ce — TF, Cg — oTT, Clo — TJiTs ; • • • l^Z.Z^ 



The g-Gaussian (12.251) contains two adjustable parameters, q and h. For the 
best fit we will adopt the criterion that the coefficients of the quadratic and 
quartic terms in the expansions (I2.25P and (12.260 coincide. This leads directly 



to the relation C4/C2 



q). Working out the remaining algebra using 



(E^SD, ^TM . and dOTj) . we obtain 



Ip 



4vr(l - Ip) 



1-p 



(2.28) 



The first one of these equations is precisely the result (12. 9p conjectured by 
Thistleton et al. 1121 on the basis of their numerical simulations! As an 



example, we show in figured] the exact solution (12.230 for p = ^ (solid line) 
together with its g-Gaussian approximation (dotted line), which according 
to (12.280 has g = — §• The difference between the two is of the order of the 
thickness of the dotted line and just barely visible to the eye. 

To get an analytic idea of the quality of the fit (I2.28P we consider the ratio 
of the g-Gaussian approximation to the true solution. /^From (12.250 - (12.281) 
we deduce that 



r{u) 



p 



2-p 



exp 



2(1 -P) 
P 



e2n7r U 



(2.29) 



n=3 
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where we set 

e2n = n-'{iy~' -C2n, (n = 1 , 2, . . .) • (2.30) 

By construction £2 = £4 = 0. It is remarkable that for the next few coeffi- 
cients the terms in the difference (12.301) largely cancel. Evaluation gives, in 
particular, 

= tI = 0.01481, 68 = ^= 0.03386, eio = ^fffg = 0.05799. (2.31) 

The smallness of these coefficients explains analytically why the g-Gaussian 
fit dOSD is so excellent. 

Finally, we observe that with the criterion employed here the two func- 
tions V and Gq do not exactly coincide for U = 0. Using the explicit expres- 
sion for Cq we find for their ratio in the origin 

^ = hhr f-^V' = 3(2 -p) [nj-DY 

V{0) '\2-p) 2(3-2p)r(A-l)r(|;) 

= l-^p2 + ..., (2.32) 

where the last line represents a small-p expansion. Again we see that the 
deviation from unity is very small. Hence our analysis confirms the high 
quahty of the fit. 

2.6 Remark on the general case (12.211) 

The general case that we set out to consider in section 12.21 is specified by an 
arbitrary marginal one- variable distribution Pi{x). The corresponding limit 
distribution 'P{U) found in fl2.2ip is then obtained in terms of a function 
!/,([/) that follows from Pi(X) via relations fIZTU]) . fimi) . fl^TSD . and flZ^ . 
Although V{U) is not a g-Gaussian in the special case where Pi{x) is a 
block function, it is likely that an appropriate choice of Pi{x) will lead to 
a Vt:{U) for which V{U) is a g-Gaussian. The conditions for this to happen 
constitute a rather complicated set of equations that we will not study here. 
The important point is that in the same way not only a g-Gaussian but any 
limit function V{U) out of a large class can be obtained; there is no indication 
that g-Gaussians play a special role. 
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3 Second example: the MTG sum 



3.1 A special probability law 

In an earlier attempt to construct a g-Gaussian distributed sum of variables, 
Moyano, Tsallis, and Gell-Mann (MTG) [13] proposed the expression 

V+n-Nl n \ P 



i=N-n 



i?^(n) = y (-1)^+"-^ . \ ^ — (3.1) 

^ \^ + n-NJ ^,_^,_l)pp]y > 



which depends on two parameters p and p [T5]. The RN{n) satisfy the nor- 
malization 



i?^(n) = l (3.2) 

n=0 ^ ^ 

and (^)i?jv(n) can be interpreted as the probability that the sum of N 
Boolean random variables, correlated in a special way (but remaining im- 
plicit), be exactly equal to n. In this second example the starting point (13.11) 
is itself clearly motivated by the formalism of g-statistical mechanics: the 
last factor in the summand represents a so-called 'g-product' [13] with p in 
the role of 1 — g; in the singular limit p — it reduces to 

Moyano et al. [13] were again interested in the large- iV limit with n scaling 
as n = yN, hence the limiting probability distribution 

n{y) = lim N(^^RM{yN), (3.3) 



N^oo \yN. 

where the first factor comes from the scaling between n and y. They 
evaluated the sum (13. ip numerically for typical values of p and p and for 
up to 1000, taking extraordinary precautions to deal with the cancellations 
caused by the alternating signs. On the basis of their numerical results 
Moyano et al. [13] conjecture that IZi^y) is a double-branched g-Gaussian, 
the two branches joining in the center and having the same g but slightly 
different widths. After fitting they arrive at the relation 

(3.4) 

1-p 

The functional dependence of (13.41) on p bears resemblance to that of g in 
the first example, equation (12. 9p . 

For simplicity we will restrict ourselves in what follows to the case p = | , 
where some rewriting turns expression (13. ip into 

1 . ^ I n 



Rj^(n) = - -> (-1)"-^ r- (3.5) 

2(l-2-P)^^ ^ \j)[N-j + {2P-l)-^f. ^ ' 

Below we will show that it is possible, again, to evaluate this sum and find 
the limit function 'R-{y) analytically. 
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3.2 The limiting law n{y) 
Upon employing in (13.51) the integral representation 

1 1 



Xf' T{k) 



daa'^-V"^ (3.6) 



withX = N-j + {2P-l)-^ and A; = 1/pwefind that the sum on j represents 
the binomial expansion of (e" — 1)". We get 

R^^n) = ^ daa^^e-^^^-i^^'^-^^^") (3.7) 
Jo 

with 

= 2r(i)(l-2-0^ 

/(«) = a-ylog(e--l). (3.8) 

The function /(a) has its maximum for a = a* = log(l — y)- In the large- 
limit the integral (13. 7p is easily evaluated by the saddle point method with 
the result 

where the dots indicate terms that vanish as ^ oo and we recall that 
y = n/N. This should now be combined with the large- iV expression of the 
binomial coefficient, 

= [2nNy{l - y)] Q-my^ogyHi-y) Mi-y)] [ i + . . . ] ^ (3.10) 

which results from straightforward application of Stirling's formula. Com- 
bining (1S3D with ([S3]) and dXTU]) and setting 

2 — 2^ 

an = (3.11) 

we obtain 

n{y) = (1 - yr^ [ - log(l -y)]"^, {0<y< 1). (3.12) 

This is our result for the distribution 7l{y) studied by Moyano et al. [13] in 
their special case p = \. Again, it is not the double-branched g-Gaussian 
that the authors [T3] conjectured it to be (and which would have required it 
to be nonanalytic at its center). 
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Figure 2: Comparison of two relationships between q and p. Solid line: equation 
(j3.15p obtained analytically in this work. Dotted line: equation (|3.4p conjectured 
by Moyano et al. [13j on the basis of numerical work. The difference is hardly 
visible to the eye. 
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Figure 3: Solid line: the exact limit function TZ{y) given by equation (13.121) for 
the case P = jq- Dotted line: the g-Gaussian approximation to this function for 
g = — I (see text). We improved the fit visually by not centering the g-Gaussian 
at the exact maximum y* but shifting it very slightly to the left. The difference 
between the two curves is small but visible. 
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3.3 Discussion of 'R'{y) 



The exponents and (1 — p)/p are both positive for < p < 1, which is 
the interval of interest. Hence TZiy) is bell-shaped and vanishes in ?/ = and 
y = 1. One may verify explicitly the normalization 7l{y)dy = 1. As had 
to be expected, TZ{y) has no point of symmetry. It takes its maximum for 
y = y*p given by 

y; = l- exp(-^). (3.13) 

^ \ pap J 

We consider again two limiting cases of interest. For p — one shows 
after some algebra that TZiy) approaches a Gaussian centered at ?/ = | and 
of variance p(log2)^. For p ^ 1 we have Ai = 1 and ai = 0, whence it 
follows that TZiy) is the block function TZ{y) = 1 on the interval < y < 1. 
When p ^ 1 the maximum yl of the distribution has the nontrivial limit 
yl = limp^i y* = 1 — exp ( — (21og2)~^) = 0.5139..., which is different from 
1 in spite of TZiy) being symmetric for p = 1. In summary, for < p < 1, 
the sum TZ{y) varies between a Gaussian and a block function distribution, 
just as in the first example. 



3.4 Best g-Gaussian fit to 7l{y) 

We will again find the g-Gaussian that best fits the exact result (13.121) . Upon 
expanding the latter around its maximum fl3.13p one finds 

oo 

n{y) = exp (- ^D,.(y - y;yy (s.u) 

which defines the coefficients Dj. The best q will be determined as a function 
of p in the same way as above, that is, from the ratio of the coefficients of 
the quadratic and the quartic terms according to |(1 — q) = D4^/D\. This 
procedure 'averages' over the asymmetry due to the cubic and higher order 
odd terms in (13.141) by simply neglecting them. After doing the algebra one 
finds [15] 

, 11(1 -p)^- 12(1 -p)pap + 6p^a^ 
which has the small-p expansion 

q = 1- [^(log2)2-41og2 + 2]p- [ll(log2)=^-f(log2)2 + 2]p2 + ... 
= 1 - 0. 98907.. .p- 1.01889.. .p2 + ... . (3.16) 

It appears from (13.151) that in this second example we do not reproduce 
the guess (13. 4p based by MTG [13] on their numerical evaluation of Tiiy). 
Nevertheless, the two relations between q and p, equations (13.41) and (I3.15p . 
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are in very close numerical agreement, as shown in figure [2l Concomitantly 
expansion (13.161) . when compared to the expansion q = 1— p— p^+. . . obtained 
from (13. 4p . shows a close numerical agreement between the coefficients. 

In figure [21 finally, we show the limit function TZ{y) for p = together 
with its g-Gaussian approximation, which according to equation (13. 4p has 
q = — |. No visible difference results if instead one were to take q = —1.36074 
as would follow from (13.151) . Hence in this second example once more, things 
conspire such that it becomes extremely difficult to distinguish the true curve 
from its g-Gaussian approximant. 

4 Final remarks 

Probability distributions shaped as g-Gaussians, if and when they occur, are 
interpreted as a manifestation of a g-generalized statistical mechanics [1] . We 
have investigated the emergence of g-Gaussians as probability distributions 
of sums of strongly correlated variables. The possible role of g-Gaussians in 
a wider context has remained outside the scope of this work. 

The literature examples discussed in this note illustrate the difficulty 
of identifying g-Gaussians. Two cases where g-Gaussian distributions had 
been thought to occur, turn out to correspond to analytically very different 
functions. The very close numerical agreement does not detract from this 
distinction of principle. Although the first example does show a way of fine- 
tuning correlations between random variables such that their sum becomes 
g-Gaussian distributed, this method works for any one out of a large class 
of distributions and does not hint at a special role for g-Gaussians. The 
probability laws constructed by this kind of fine-tuning are limit distributions 
in the plain sense that they appear when the number N of terms in the sum 
tends to infinity. They are not, however, or at least have not been shown 
to be, attractors, and this raises the question of why such finely adjusted 
correlations would occur in physics. Efforts to let g-Gaussians play the role 
of attractors within the framework of a g-generalized central limit theorem 
[T6] have not yet led to examples that we can subject to analysis. 

The indeniably interesting properties of g-Gaussians make it worthwile to 
investigate the role they play in many-body systems in statistical mechanics. 
However, from the examples studied in this work we draw the conclusion 
that, in the absence of convincing theoretical arguments, one should ex- 
ercise extreme caution when interpretating non-Gaussian data in terms of 
g-Gaussians. 
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